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ABSTRACT 

The contribution of this paper is in suggesting an 
analysis and design of a control system with variable 
parameters. By applying the recommended by the 
author method of the Advanced D-partitioning the 
system’s stability can be analyzed in details. The 
method defines regions of stability in the space of the 
system’s parameters. The designed controller is 
enforcing desired system performance. The suggested 
technique for analysis and design is essential and 
beneficial for the further development of control 
theory in this area. 
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I. INTRODUCTION 

Control systems performance must be insensitive to 
parameter variations. In the process of design such a 
control system, it is important to determine the 
regions stability, related to the variation of the system 
parameters. The suggested by the author method, 
dealing with the effects of parameters variations on 
the system’s stability, is classified as Advanced D- 
partitioning [1], [2], [3]. It is an efficient tool for 
system stability analysis in case of variation of any of 
the system’s parameters. It is further upgraded in this 
paper and can be used for simultaneously varying 
parameters. This research is also suggesting a method 
for design of a controller, by applying forward-series 
compensation. It can suppress the influence of any 
parameters variations of the control system. 
Innovation is demonstrated in the unique property of 
the designed controller that can operate effectively for 
variations of any one of the system’s parameters 
within prescribed limits. The design of the controller 
is based on the ITAE criterion [4], [5], For higher 
order systems, a pair of dominant poles represents the 
system dynamics. The relative damping ratio if of the 


system is taken as a performance objective for the 
optimization design. 


II. System with Simultaneously Variable 
Parameters 

The Advanced D-Partitioning analysis in case of two 
simultaneously variable parameters [2], [4], [6] can be 
demonstrated for a control system of the armature- 
controlled dc motor and a type-driving mechanism. 
The gain and one of the time-constants are uncertain 
and variable. The open-loop transfer function of the 
system can be presented as: 


G P0 ( s ) - 


_ K _ 

(l + ?y)(l + 0.5s)(l + 0.8s) 


( 1 ) 


The characteristic equation of the unity feedback 
system is derived as: 

K + (1 + 7s)(l + 0.5s)(l + 0.8s) = 0 (2) 


By substituting s = jco equation (2) is modified to: 

K = -1 + (1 3T + 0.4)<y 2 + ja)(0ATa) 2 -13-T) (3) 


Since the gain may have only real values, the 
imaginary term of equation (3) is set to zero. Then: 

(4) 

OAT v 7 


The result of (4) is substituted into the real part of 
equation (3), from where: 


1.3r»+14jW+a52 . 3Jsr+4ias +L3 (5) 

OAT T 


The D-Partitioning curve K = f (T) is plotted with the 
aid of the following code: 

» T = 0:0.1:5; 

» K = 3.25.*T+4.225+1.3./T 
K = 
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Columns 1 through 10 

Inf 17.5500 11.3750 9.5333 8.7750 8.4500 

8.3417 8.3571 8.4500 8.5944 

Columns 11 through 20 

8.7750 8.9818 9.2083 9.4500 9.7036 9.9667 
10.2375 10.5147 10.7972 11.0842 


Columns 21 through 30 

11.3750 11.6690 11.9659 12.2652 12.5667 

12.8700 13.1750 13.4815 13.7893 14.0983 


Columns 31 through 40 

14.4083 14.7194 15.0313 15.3439 15.6574 

15.9714 16.2861 16.6014 16.9171 17.2333 

Columns 41 through 50 

17.5500 17.8671 18.1845 18.5023 18.8205 

19.1389 19.4576 19.7766 20.0958 20.4153 

Column 51 
20.7350 
» Plot (T, K) 


D-Partitioning by Two Variable Parameters (System Type 0( 



Fig. 1. Advanced D-Partitioning in terms of two 
variable parameters 


The D-Partitioning curve K = f (T) defines the 
border between the region of stability D (0) and 
instability D(l) for case of simultaneous variation of 
the two system parameters. Each point of the D- 
Partitioning curve represents the marginal values of 
the two simultaneously variable parameters. This is a 
unique advancement and an innovation in the theory 
of control systems stability analysis. The 
demonstration of the system performance in case of 
variation of the time-constant T is done at gain set to 
K = 10. When 0 < T < 0.25 sec and T > 1.5 sec the 


system is stable. But it becomes unstable in the range 
0.25 sec < T < 1.5 sec. It is also obvious that the 
system performance and stability depends on the 
interaction between the two simultaneously varying 
parameters. If K < 8.3417, the system is stable for 
any value of the T. Higher values of K (K = 12, K = 
14), enlarge the range of T at which the system will 
fall into instability. 


III. Design of a Controller for Systems with 
Simultaneously Variable Parameters 

The open-loop transfer function of the plant Gpo(s) is 
modified and now presented in equation (6) 
considering the two variable parameters that are the 
system’s gain K and time-constant T. Initially, it is 
suggested that the gain is set to K = 10, while T is 
variable [5], [7], [8], 

K 


G P0 (s) = 


(1 + 73)(1 + 0.55)(1 + O. 85 ) 
_ K _ 

0.47s 3 + (1 3T + 0.4)s 2 +(1.3 + T)s +1 


( 6 ) 


The robust controller consists of a series stage Gsoi s) 
and a forward stage Gfo( s). An integrating stage G/o( s) 
is also included in the controller as seen from Figure 
2 . 



Fig.2. Robust controller incorporated into the control 

system 

Initially, the plant transfer function Gpo (s), as a 
standalone block, is involved in a unity feedback 
system having a closed-loop transfer function 
presented as: 


G cl (s) = 


_ K _ 

(l + Ts)(l + 0.5s)(l + 0.Ss) + K 


K 


(7) 


OATs 3 +(1.3r + 0.4> 2 +(1.3 + T)s + 1 + K 


Equation (7) is used as a base in the design strategy 
for constructing the series stage of the robust 
controller. It has the task to place its two zeros near 
the desired dominant closed-loop poles, that satisfy 
the condition C = 0.707 [9], [10], These zeros will 
become the dominant poles of the unity negative 
feedback system, involving the cascade connection of 
the series controller stage Gso(s), the integrator Gio(s) 
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and the plant Gpo(s). If initially the gain is set to K = 
10, the optimal value of the time-constant T, 
corresponding to the relative damping ratio C = 0,707 
of the closed-loop system, is determined by the code: 

» T= [20:0.01:35]; 

» For n=l: length (T) 

G_array(:,:,n)=tf([10],[0.4*T(n)(1.3*T(n)+0.4) 
(1.3+T(n)) ll]);end 
» [y, z] = damp (G_array); 

» [y, z] = damp (G_array); 

» Plot (T,z(l,:)) 

Relative Damping Ratio in Terms of Time-Constant 



As seen from Figure 3, the relative damping ratio is 
C = 0.707 when the time-constant is T = 27.06 sec. 

By substituting the values T = 27.06 sec and K = 10 in 
equation (7), the transfer function of the closed-loop 
system becomes: 

10 


G cl {s)~ 


10.824s 3 + 35.578s 2 + 28.36s +11 


( 8 ) 


The series robust controller zeros can be placed at the 
approximated values -0.5 ± j0.5. Therefore, the 
transfer function of the series robust controller Gso(s) 
is: 

(s + 0.5 + ;'0.5)(s + 0.5 - ;0.5) 


G. n (s) = 


0.5 


s + s + 0.5 
0.5 


(9) 


An integrating stage G/o(s) is added to eliminate the 
steady-state error of the system. It is connected in 
cascade with the series controller. Then, the transfer 
function of the compensated open-loop control system 
will be as follows: 


GqlG) — G / fl(s)C7< ;0 (.s , )C7p 0 (s) — 


'/o 


so' 


’P o' 


K(s z + s + 0.5) 


( 10 ) 


0.5s(l + 7s)(l + 0.5s)(l + 0.8s) 


When Gol( s) is involved in a unity feedback, its 
closed-loop transfer function is determined as: 

iv ,> 

,(4rrc 


K(s 2 + s + 0.5) 


CL 


0.5s(l + 7s)(l + 0.5s)(l + 0.8s) + 
+ K(s 2 + s + 0.5) 


(ID 


It is seen from the equation (11) that the closed-loop 
zeros will attempt to cancel the closed loop poles of 
the system, being in their area. This problem can be 
avoided if a forward controller Gfo( s) is added to the 
closed-loop system, as shown in Figure 2. The poles 
of Gfo(s) are designed to cancel the zeros of the 
closed-loop transfer function Gcl(s), as shown in 
equation (12): 

0.5 


G fo ( s ) 


f v ■ 0.5 


( 12 ) 


The assessment of the system proves that the relative 
damping ratio becomes C = 0.707, when the time- 
constant is T = 27.06 sec, resulting in system’s desired 
closed-loop poles -0.466 ± j0.466. These outcomes 
are determined from the code: 

» GCL0=tf([10],[10.824 35.578 28.3611]) 

» damp (GCL0) 


Finally, the transfer function of the total compensated 
system is derived considering the diagram in Figure 2. 


G T q (s) — G fo G clos (s) — 
K 


0.2 Ts 4 + (0.65T + 0.2)s 3 + 


+ (0.65 + 0.57> 2 + 0.5j + K(s 2 + s + 0.5) 


(13) 


Eigenvalue 

-4.64e-001 + 4.64e-001i 
-4.64e-001 - 4.64e-001i 
-2.36e+000 


Damping Freq. (rad/s) 
7.07e-001 6.56e-001 

7.07e-001 6.56e-001 

1.00e+000 2.36e+000 


IY. Performance of the Compensated System 

The system is tested for insensitivity to variations of 
its gain K and its time-constant T. Comparison of its 
performance is done before and after applying the 
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robust compensation. Initially at system gain is set to 
K = 10 and three different values of the time-constants 
are set successively to T= 0.1 sec, T = 0.8 sec and T = 
2 sec and are substituted in equation (6). 

The case of T = 0.8 sec, corresponds to the region 
D(2) and definitely to an unstable control system. The 
cases of T =0.1 sec and T = 2 sec, reflect regions 
D1(0) and D2(0) accordingly and are related to a 
stable control system. 

The transient responses of the system before applying 
the robust compensation are illustrated in Figure 4 and 
Figure 5 and are achieved by the following codes: 

» Gp001= tf([010],[0.04 0.53 1.41]) 

» Gp02= tf([010],[0.8 3 3.3 1]) 

» Gp0fb001=feedback (GpOOl, 1) 

» Gp0fb02=feedback (Gp02,1) 

» Step (GpOfbOOl, Gp0fb02) 


Step Response of the Original System Type 0 
(Cases T = 0.1 sec. T = 2.0 sec| 



Fig.4. Step responses of the original control system 
(T = 0.1 sec, T = 2sec at K = 10) 

» Gp008=tf([010],[0.32 1.442.11]) 

» Gp0fb008=feedback (GpOOl,1) 

» step (Gp0fb008) 


Step Response of the Original System Type 0 
(Case T = 0.8sec) 



Fig.5. Step responses of the original control system 
(T = 0.8sec, at K = 10) 


The compensated system is also examined for 
robustness in the time-domain. Considering the same 
values of the time-constant, as those used for the 
assessment of the original system, T= 0.1 sec, T = 0.8 
sec and T = 2 sec at system gain K = 10 and 
substituting them in equation (13), the following 
transfer functions are obtained: 


G t As) 


10 


T = 0.1 


0.04s 4 + 0.53s 3 + 21,4s 2 + 21s +10 


(14) 


^Yo ( S ^ T 


0.8 


10 

0.32s 4 + 1.44s 3 + 22.1s 2 + 21s +10 


(15) 


G T0 (s) = - (16) 

T ~ 2 0.8s 4 + 3s 3 + 23.3s 2 + 21s+ 10 

The step responses, representing the time-constant 
variation of the robust system, shown in Figure 6, are 
determined by the code: 

» GT01=tf([10], [0.04 0.53 21.4 2110]) 

» GT08=tf([10], [0.32 1.44 22.1 2110]) 

» GT20=tf([10], [0.8 3 23.3 2110]) 

» step (GT01, GT08, GT20) 


Step Response of the System Type 0 with a Robust Controller 



Fig.6. Step responses of the system with a robust 
controller (T = 0. 1 sec, T = 0. 8 sec, T = 2 sec at K = 

10 ) 


As seen from Figure 6, due to the applied robust 
controller, the control system becomes quite 
insensitive to variation of the time-constant T. The 
step responses for T = 0. 1 sec, T = 0. 8 sec and T = 
2 sec coincide. 


Since the system is with two variable parameters, now 
the gain will be changed, applying: K = 5, K = 10, K = 
20, while keeping the system’s time-constant at T = 
0.8 sec. For comparison of the system’s insensitivity 
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to the gain variation before and after applying the 
robust compensation, initially the suggested values as 
shown above are substituted in equation (6). The 
transient responses of the original system are 
illustrated in Figure 7 and 8 and are achieved by the 
following codes: 

» Gp05=tf([0 5],[032 1.44 2.11]) 

» Gp010=tf([0 10],[0.321.442.11[) 

» Gp0fb5=feedback (Gp05,1) 

» GpOfblO=feedback (GpOlO, 1) 

» Step (Gp0fb5, GpOfblO) 

» Step (Gp0fb5, GpOfblO) 


Step Response of the Original System Type 0 
|Cases K = 5. K = 10) 



Fig.7. Step responses of the original control system 
(K = 5, K = 10 at T = 0.8sec) 


» Gp020=tf([0 20],[0.32 1.44 2.11[) 
» Gp0fb5=feedback (Gp05,1) 

» GpOfblO=feedback (GpOlO, 1) 

» Gp0fb20=feedback (Gp020,1) 

» step (Gp0fb20) 


Step Response of the Original System Type 0 
(Case K = 20) 



Fig. 8. Step response of the original control system 
(K = 20 at T = 0.8sec) 


It is obvious that the cases of K = 10 and K = 20, 
correspond to an unstable original control system. 

Next, the variable gain K = 5, K = 10, K = 20 will be 
applying to the robust compensated system, keeping 
the system’s time-constant at T = 0.8 sec. These 
values are substituted in equation (13). As a result, the 
following outcomes are delivered: 

5 

G ro ( 5 ) =- (17) 

K ~ 5 0.32s 4 + 1.44s 3 + 12.1s 2 + 11s + 5 

G T0 (s) =--- (18) 

K ~ W 0.32s 4 + 1.44s 3 + 22.1s 2 + 21s+ 10 

G ro (s) =--- (19) 

K ~ 20 0.32s 4 + 1.44s 3 + 42.1s 2 + 41s+ 20 

To compare the system robustness before and after the 
robust compensation, the step responses for the three 
different cases, representing the gain variation of the 
robust system, are plotted in Figure 9 with the aid of 
the code as shown below: 

» GTK5=tf([2.5],[0.16 0.72 6.05 5.5 2.5]) 

» GTK10=tf([5],[0.16 0.72 11.05 10.5 5[) 

» GTK20=tf([10],[0.16 0.72 21.05 20.510]) 

» step (GTK5, GTK10, GTK20) 


Step Response of the System Type 0 with a Robust Controller 



Figure9. System’s step responses with the robust 
controller (K = 5, K = 10, K = 20 at T = 0.8 sec) 


Again, the system became robust. An additional zero 
included into the series robust controller stage can 
further improve the rise time of the system’s step 
response. 

An average case is chosen with K = 10 and T = 0.8 
sec for the assessment of the system’s performance 
after the application of the robust controller. This case 
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will differ insignificantly from the other cases of the 
discussed variable K and T. The performance 
evaluation is achieved by following code: 


» GT10=tf ([10], [0.32 1.4422.1 2110]) 
» damp (GT10) 


Eigen value Damping 

-4.91e-001 + 4.89e-001i 7.09e-001 
-4.91e-001 - 4.89e-001i 7.09e-001 
-1.76e+000 + 7.88e+000i 2.18e-001 
-1.76e+000 - 7.88e+000i 2.18e-001 


Freq. (rad/s) 
6.93e-001 
6.93e-001 
8.07e+000 
8.07e+000 


It is seen that the relative damping ratio enforced by 
the system’s dominant poles is C = 0.709, being very 
close to the objective value of C = 0.707. This 
insignificant difference is due to the rounding of the 
desired system’s poles to -0.5 ± j0.5, during the 
design of the series robust controller stage. 

Y. Conclusions 

The D-Partitioning analysis is further advanced for 
systems with multivariable parameters [8], [9], [10], 
The Advanced D-partitioning in case two variable 
parameters is demonstrating the strong interaction 
between the variable parameters. Each point of the D- 
Partitioning curve represents the marginal values of 
the two simultaneously variable parameters, being a 
unique advancement and an innovation in the theory 
of control systems stability analysis. 

The design strategy of a robust controller for linear 
control systems proves that by implementing desired 
dominant system poles, the controller enforces the 
required relative damping ratio and system 
performance. For systems Type 0, an additional 
integrating stage ensures a steady-state error equal to 
zero. 

The designed robust controller brings the system to a 
state of insensitivity to the variation of its parameters 
within specific limits of the parameter variations. The 
experiments with variation of different parameters 
show only insignificant difference in performance for 
the different system conditions [10], [11], [12], 

For the discussed case, the system becomes quite 
insensitive to variation of the time constant within the 
limits 0.1 T < T < 10 T. The system is quite insensitive 
to variations of the gain K within the limits 0.5 K < K 
< 5 K. Insignificant step response difference is 
observed also if the experiment is repeated with 


different variation of the gain and different variation 
time-constants values. 

Since the design of the robust controller is based on 
the desired system performance in terms of relative 
damping, its contribution and its unique property is 
that it can operate effectively for any of the system’s 
parameter variations or simultaneous variation of a 
number of parameters. This property is demonstrated 
by the comparison of the system’s performance before 
and after the application of the robust controller. Tests 
demonstrate that the system performance in terms of 
damping, stability and time response remains robust 
and insensitive in case of any simultaneous variations 
of the gain and the time-constant within specific 
limits. The suggested analysis and design is beneficial 
for further advancement of control theory in this field. 
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